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We consider scalar-tensor theories of gravity in an accelerating universe. 
The equations for the background evolution and the perturbations are given 
in full generality for any parametrization of the Lagrangian, and we stress 
that apparent singularities are sometimes artifacts of a pathological choice of 
variables. Adopting a phenomenological viewpoint, i.e., from the observations 
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back to the theory, we show that the knowledge of the luminosity distance as 
a function of redshift up to z ~ (1 — 2), which is expected in the near future, 
severely constrains the viable subclasses of scalar-tensor theories. This is due 
$h ' to the requirement of positive energy for both the graviton and the scalar 

partner. Assuming a particular form for the Hubble diagram, consistent with 
present experimental data, we reconstruct the microscopic Lagrangian for 
various scalar-tensor models, and find that the most natural ones are obtained 
if the universe is (marginally) closed. 

PACS numbers: 98.80. Cq, 04.50.+h 
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I. INTRODUCTION 

Recently, there has been a lot of interest in cosmological solutions in the presence of a 
cosmo logical constant, when the latter is significant compared to the present total energy 
density of the universe. Indeed, the Hubble diagram based on observations of type la su- 
pernovae up to a redshift z ~ 1 seems to imply that our universe is presently accelerating 
|l]||. These data, when combined with the observed location of the first acoustic peak of 
the CMB temperature fluctuations, favor a spatially flat universe whose energy density is 
dominated by a "cosmological constant" -like term. The flatness of the Universe is corrobo- 
rated by the latest Boomerang and Maxima data [§|f|], in accordance with the inflationary 
paradigm, though a marginally closed Universe is still allowed by the position of the first 
acoustic (Doppler) peak at I ~ 200. A significant cosmological constant may help in re- 
solving the dark matter problem - for dustlike matter alone observations seem to imply 
Q m ~ 0.3 - and in reconciling flat Cold Dark Matter (CDM) models with observations in 
the framework of ACDM models. Finally, a cosmological constant is an elegant way to allow 
a high Hubble constant H with h = H /(100 km s" 1 Mpc -1 ) « 0.65 and a sufficiently old 
universe to > HGyr || (see also, e.g., || for a recent comprehensive review and references 
therein). 

Therefore, this interpretation, if confirmed by future observations, constitutes a fun- 
damental progress towards the solution of the dark matter problem and the formation of 
large-scale structure in the Universe out of primordial fluctuations generated by some infla- 
tionary model. That is certainly what makes it so appealing and gives it, maybe somehow 
prematurely, the status of new paradigm. A striking consequence for our Universe is then 
its present acceleration, for a large range of equations of state J7|. 

Of course, from the point of view of particle physics, a pure cosmological constant of 
the order of magnitude A w 3 x 10 -122 c 3 /(HG), interpreted as the vacuum energy, is ex- 
tremely problematic. This is why attempts were made to find some alternative explanation 
to the origin of the acceleration under the form of some scalar field $ (sometimes called 
quintessence ||, "A"-field, etc.) whose slowly varying energy density would mimic an ef- 
fective cosmological constant. This is very reminiscent of the mechanism producing the 
inflationary phase itself with the fundamental difference that this scalar field, which does 
not have to be a priori the inflaton, is accelerating the expansion today, therefore at a much 
lower energy scale. This of course has problems of its own as this effective cosmological 
constant term started dominating the universe expansion only in the very recent past (the 
so-called "cosmic coincidence" problem). Indeed, the energy density of the field $ must re- 
main subdominant at very early stages and come to dominate in the recent past only. Hence, 
specific evolution properties are required to meet these constraints and were indeed shown 
to hold for particular potentials, partly alleviating the problem of the initial conditions. For 
inverse power-law potentials the energy density of the scalar field was shown to decrease less 
rapidly than the background energy density so that it can be negligible in the early universe 
and still come to dominate in the recent past M. For exponential potentials, pH|J5|] the scalar 
field energy density has the very interesting behavior that it tends to a fixed fraction of the 
total energy density, these are the so-called "tracker solutions" . Hence a pure exponential 
potential is excluded if data confirm that the energy density of the scalar field is dominating 
today, as this fraction had to be small at the time of nucleosynthesis. A slightly different 



potential is proposed in [I 1 1 and a classification of the scaling behavior of the scalar field 



for various potentials has been given in ||12|| . Hence, though a minimally coupled scalar field 
is an attractive possibility, some degree of fine tuning still remains in the parameters of the 
potential [0,0]. 



If one admits that it is some minimally coupled scalar field which plays the role of an 
effective cosmological constant while gravity is described by general relativity, the question 
immediately arises: What is the "right" potential £/($) of this scalar field? In a recent work 
by Starobinsky |I4] , the following "phenomenological" point of view was adopted: Instead 



of looking for more or less well motivated models, like the interesting possibilities discussed 
above, it is perhaps more desirable to extract as much information as possible from the 
observations (a similar approach can also be adopted to reconstruct the inflaton potential) 
in order to reconstruct the scalar field potential, if the latter exists at all. Cosmological 
observations could then be used to constrain the particle physics model in which this scalar 
field is supposed to originate. In the context of general relativity plus a minimally coupled 
scalar field, it was shown that the reconstruction of U(<&) can be implemented once the 
quantity Dl(z), the luminosity distance as a function of redshift.is extracted from the ob- 
servations [|T4| , p~5|| , something that is expected in the near future.tl The SNAP (Supernovae 
Acceleration Probe) satellite will notably make measurements with an accuracy at the per- 
cent level up to z ~ 1.7. Of course, in this way only the recent past of our Universe, up to 
redshifts z ~ (1 — 2) (for reference, we will push some of our simulations up to z ~ 5), is 
probed and so the reconstruction is made only for the corresponding part of the potential. 
Crucial information is therefore gained on the microscopic Lagrangian of the theory through 
relatively "low" redshift cosmological observations. 

A further step is to generalize the same mechanism in the framework of scalar-tensor 
theories of gravity, sometimes called "generalized quintessence" . The usual minimally cou- 
pled models are certainly ruled out if, for example, it turns out that this component of the 
energy density obeys an equation of state p = wp with w < — 1 (p > 0). Strangely enough, 
such an unexpected equation of state which in itself implies new physics, is in fair agreement 
with the observations |17| . Also the inequality dH 2 (z)/dz > 3fi mi o-ffo(l + z) 2 must hold for 
a minimally coupled scalar field, hence its violation would force us to consider more compli- 
cated theories, possibly scalar-tensor theories. There are also strong theoretical motivations. 
These theories, in which the scalar field participates in the gravitational interaction, are the 
most natural alternatives to general relativity (GR). Indeed, scalar partners to the graviton 
generically arise in theoretical attempts at quantizing gravity or at unifying it with other 
interactions. For instance, in superstrings theory, a dilaton is already present in the super- 
multiplet of the 10-dimensional graviton, and several other scalar fields (called the moduli) 
also appear when performing a Kaluza-Klein dimensional reduction to our usual spacetime. 
Moreover, contrary to other alternative theories of gravity, scalar-tensor theories respect 
most of GR's symmetries: conservation laws, constancy of (non-gravitational) constants, 
local Lorentz invariance (even if a subsystem is influenced by external masses), and they 



1 Actually, it is shown in Ref. [16] that the potential U(&) can already be reconstructed from 



present experimental data, although not yet very accurately. 



also have the capability of satisfying the weak equivalence principle (universality of free fall 
of laboratory-size objects) even for a strictly massless scalar field. Nevertheless, they can 
describe many possible deviations from GR, and their predictions have been thoroughly stud- 
ied in various situations: solar-system experiments |18|JT9|J20[| , binary-pulsar tests ||T8| . pT)| , |2l| . 
gravitational- wave detection f22| . |23|j . Finally these scalar-tensor theories could play a crucial 



role in the very early universe, for example in the Pre Big Bang inflationary model (see e.g. 

0). ' 

Thus, in this work we are investigating the possibility to have an accelerating universe 
in the context of scalar-tensor theories of gravity instead of pure GR. This has indeed 
attracted a lot of interest recently and such cosmological models have been studied and 
possibly confronted with observations like CMB anisotropies, or the growth of energy density 
perturbations (see for instance [p5| , p6| , p7| , p8| , p9| , p0| , pT1 , p2| , p3| , p4| , p5| ) . However, we emphasize 



once more that the central point of view adopted here, in analogy with Starobinsky ||14|| , 
is to constrain the model with the experimental knowledge of the Hubble diagram up to 
z ~ (1 — 2). This is precisely why use of the redshift z as basic variable is crucial for our 
purpose: Quantities like H (z) are directly observable, in contrast to, say,H H(t) or H(Q). For 
instance, we have access to H (z) through the direct measurement of the luminosity distance 



in function of redshift Dl(z). In a recent letter [3B|, it was shown that the knowledge of 
both H(z) and S m (z) is sufficient to reconstruct the full theory (again, in the range probed 
by the data). This means that we do not choose any specific theory a priori, but instead we 
reconstruct whatever theory possibly realized in Nature. 

As we will see, the knowledge of H(z) on its own, though insufficient in order to fully 
reconstruct a scalar-tensor theory unless one makes additional assumptions, turns out to 
be already very constraining when subclasses of models are considered. This is particularly 
interesting because it means that cosmological observations at low redshifts implying an 
accelerated expansion might well give new constraints on scalar-tensor theories. We will 
show that this is indeed the case. 

Throughout the paper, we use natural units for which h — c — 1, and the signature 
(-+++), together with the sign conventions of [|37]]. In Section II, we introduce the general 
formalism of scalar-tensor theories of gravity and their different parametrizations. In Section 
III, we briefly review the severe experimental restrictions imposed on these theories today. 
In Section IV, we consider FRW universes in the framework of scalar-tensor gravity and 
we give the equations for the different parameterizations. In Section V, we review the full 
reconstruction problem. In Section VI, we give a detailed study of subclasses of models, 
which are investigated using the background equations. Finally, in Section VII, our results 
are summarized and discussed. 



2 The function H{t) can be obtained from the knowledge of H(z) thanks to the relation t 
- J dz/[(l + z)H(z)], but the directly observable quantity is H{z). 



II. SCALAR-TENSOR THEORIES OF GRAVITY 

We are interested in a universe where gravity is described by a scalar-tensor theory, and 



we consider the action |38] 



S = j^- J rf 4 xv^(F($) R - Z{$) g^d^d^ - 2U{$j) + S m [ip m , g^\ . (2.1) 

Here, G> denotes the bare gravitational coupling constant (which differs from the measured 
one, see Eq. (|3.5|) below), R is the scalar curvature of g^, and g its determinant. In Ref. p6| . 
we used different conventions (corresponding to the choice 87rG> = 1 in the above action); 
here, the quantity -F($) is dimensionless. This factor -F($) needs to be positive for the 
gravitons to carry positive energy. The action of matter S m is a functional of some matter 
fields ipm and of the metric g^, but it does not involve the scalar field $. This ensures that 
the weak equivalence principle is exactly satisfied. 

The dynamics of the real scalar field $ depends a priori on three functions: F(<&), Z($), 
and the potential U(Q). However, one can always simplify Z(Q) by a redefinition of the 
scalar field, so that -F($) and Z(Q) can be reduced to only one unknown function. Two 
natural parametrizations are used in the literature: (i) the Brans-Dicke one, corresponding 
to -F($) = $ and Z(&) = w($)/<l>; and (ii) the simple choice Z($) = 1 and -F($) arbitrary. 
This second parametrization is however sometimes pathological. [The derivatives of $ can 
become imaginary in perfectly regular situations; see the discussion about Eq. ( |5.6a| ) below.] 
In the following, we will write the field equations in terms of the two functions F(§) and 
Z($), so that any particular choice can be recovered easily. 

The variation of action (|2.1|) gives straightforwardly 



F($) [r^ - ^R\ = SvrG.TV + Z($) (df&dv® - ^g^(d^) 2 

+V fl d v F($) - g^DF{$) - gta ,U($) , (2.2a) 

2Z(») M = -^ B -^ («.•)» + 2§, (2.2b) 

V„T„" = , (2.2c) 

where T = T^ is the trace of the matter energy-momentum tensor T^ u = (2/\f—g) x 
SSm/Sg^. The scalar-field equation ( |2.2b| ) can of course be rewritten differently if one uses 
the trace of Eq. (|2.2aj) to replace the curvature scalar R by its source, and one gets the 



Brans-Dicke-like equation 

where 1w = 2ZF + 3(dF/d§) 2 . [In the Brans-Dicke representation where F = $ and Z = 
cij($)/$, this factor 2m reduces to the well-known expression 2cj($) + 3.] In the following, 
we will however use the form ( |2.2b| ), which will simplify considerably our calculations. 

The above equations are written in the so-called Jordan frame (JF). Since in action 
( [2.1D , matter is universally coupled to g^ u , this "Jordan metric" defines the lengths and 
times actually measured by laboratory rods and clocks (which are made of matter). All 



experimental data will thus have their usual interpretation in this frame. In particular, the 
observed Hubble parameter H and the measured redshifts z of distant objects are Jordan- 
frame quantities. 

However, it is usually much clearer to analyze the equations and the mathematical con- 
sistency of the solutions in the so-called Einstein frame (EF), defined by diagonalizing the 
kinetic terms of the graviton and the scalar field. This is achieved thanks to a conformal 
transformation of the metric and a redefinition of the scalar field. Let us call g* and ip the 
new variables, and define 

= F{$) 9fiu , (2.4a) 

= 3/ rflnF(g) \ a m_ 

" 4 V d$ ) + 2F($) [ Dj 

A{<p) = F- 1 ' 2 ^) , (2.4c) 

2V(<p) = U{$) F~ 2 ($) . (2.4d) 

Action (E. II) then takes the form 




S= ^Gj ^^^ (T " \ g * Ud ^ 9 ^ ~ V &)) + SmWml A2 &) 9U , (2-5) 

where g* is the determinant of g* , g% u its inverse, and R* its scalar curvature. Note that the 
first term looks like the action of general relativity, but that matter is now explicitly coupled 
to the scalar field tp through the conformal factor A 2 ((p). Quantities referring to the Einstein 
frame will always have an asterisk (either in superscript or in subscript), e.g. V* and □* 
for the covariant derivative and the d'Alembertian with respect to the Einstein metric. The 
indices of Einstein-frame tensors will also be lowered and raised with the Einstein metric 
g* and its inverse g% v . The field equations deriving from action (|2.5| ) take the simple form 

R% - \R*g% = SttG.t;, + 2d&d vV - g%(g?d a <pd0<p) - 2V( v ) g ; u , (2.6a) 

□V = -ArtGMv) T* + dV(<p)/d<p , (2.6b) 

VJT& = a{<p) T*d v¥ > , (2.6c) 

where 

. . din A , 

°kP) = -7— (2-7) 

dip 

is the coupling strength of the scalar field to matter sources [pJ|] , and T* = g^ u T^ u is the trace 
of the matter energy-momentum tensor T>f u = (2/y/—g*) 5S m /5g*^ v in Einstein-frame units. 
From its definition, one can deduce the relation T* = A 2 (ip) T^ with its Jordan-frame 
counterpart. 

Let us underline that the Cauchy problem is well posed in the Einstein frame [|HJ, because 
all the second-order derivatives of the fields are separated in the left-hand sides of Eqs. ( p.6[) , 
whereas they are mixed in the JF equations ( |2.2j ). Action fl2.5| ) also shows that the helicity-2 
degree of freedom is described by the fluctuations of the Einstein metric g* (whose kinetic 

6 



term is the standard Einstein-Hilbert one), and that the EF scalar ip is the true helicity-0 
degree of freedom of the theory (since its kinetic term has the standard form). On the other 
hand, the fluctuations of the Jordan metric g^ actually describe a mixing of helicity-2 and 
helicity-0 excitations, and the JF scalar $ is related to the helicity-0 degree of freedom via 
the complicated relation ( |2.4b| ), because its kinetic term in action ( |2.1| ) comes not only from 
the naive contribution Z($) (d^) 2 but also from the cross term F(Q) R. In conclusion, 
the theory can be mathematically well defined only if it is possible to write the EF action 
( |2.5| ), notably with its negative sign for the scalar-field kinetic term (so that ip carries positive 
energy). If it happens that the transformation ( |2.4| ) is singular for particular values of $, the 
consistency of the theory should be analyzed in the EF. Some singularities may be artifacts 
of the parametrization which is chosen to write action (|2.1| ), and may not have any physical 
significance. On the other hand, Jordan-frame quantities may look sometimes regular while 
there is an actual singularity in the Einstein frame (a typical example is provided when -F($) 
vanishes). In this case, the solution should be considered as mathematically inconsistent. 
In the following, we will see that the JF is better suited than the EF for our cosmological 
study, but we will always check the consistency of our results by finally translating them in 
terms of Einstein-frame quantities. 

III. KNOWN EXPERIMENTAL CONSTRAINTS 

The predictions of general relativity in weak- field conditions, and at present, are con- 
firmed by solar-system experiments at the 0.04% level PP|J4H|] . One should therefore verify 
that the scalar-tensor models we are considering are presently close enough to Einstein's 
theory. 

If the scalar field is very massive (say, if d 2 V jdp 2 is large with respect to the inverse 
of the astronomical unit), its influence is exponentially small in solar-system experiments, 
even if it is strongly coupled to matter. This situation corresponds to the particular scalar- 
tensor model considered in Ref. |4l| (namely F(&) = $ and Z(Q) = in action ( |2.i|) , but 
assuming a large enough value for d 2 U/d& 2 ). Although this situation is phenomenologically 
acceptable, it remains somewhat problematic from a field theoretical viewpoint, since the 
massive scalar would a priori desintegrate into lighter (matter) particles. 

On the contrary, if the scalar mass is small with respect to the inverse solar-system dis- 
tances, it must be presently very weakly coupled to matter for the theory to be consistent 
with experimental data. At the first post-Newtonian order (1/c 2 with respect to the New- 
tonian interaction), the deviations from general relativity can be parametrized by two real 
numbers, that Eddington jE|] denoted as {(3 — 1) and (7 — 1). In the present framework, 
they take the form [0,0,0] 



n,-1 -9 ^ {dF/d*)* m^ 

7 i_ l + a 2 ~ ZF + 2(dF/d<$>) 2 ' {6 aj 

1 a 2 da I F (dF/d$) d-y 

^~ ~ 2 (1 + a 2 ) 2 dp~~4 2ZF + 3{dF/d$) 2 rf$ ' ( ' ^ 



where the first expressions are given in terms of the Einstein-frame notation ( |2.5|) -(p77), 



whereas the last ones correspond to the Jordan-frame general representation (271). To 



simplify, the second expression of Eq. ( |3.1rj| ) has been written in terms of the derivative of 
( p. la| ) with respect to $. 



Using the upper bounds on (7 — 1) from solar-system measurements p5 |, we thus get 
the constraint 

2a 2 « (ZF) l (dF/d<S>)l < 4 x 10" 4 , (3.2) 

where an index means the present value of the corresponding quantity. On the other hand, 
the experimental bounds on {(3 — 1) cannot be used to constrain the derivative (da/d<p)o 
appearing in Eq. ( |3.1b| ), since it is multiplied by a factor a 2 , consistent with 0. Because of 



nonperturbative strong-field effects, binary-pulsar tests are however directly sensitive to this 
derivative, i.e., to the ratio — 4(/3 — l)/(7 — 1). In a generic class of scalar-tensor models, 



Refs. |2T|j2^] have obtained the bound 

(da/d(p) > -4.5 . (3.3) 



From action Eq. (|2.1| ), one can naively define Newton's gravitational constant as the 
inverse factor of the curvature scalar R: 

G N = G*A 2 = G*/F . (3.4) 

However, Gn does not have the same physical meaning as Newton's gravitational constant in 
GR. Indeed, the actual Newtonian force measured (in Cavendish-type experiments) between 
two close test masses mi and iri2 is of the form G e smim2/r 2 , where the effective gravitational 
constant reads 



C ~C A 2 (l I n 2 )- G*( 2ZF + 4(dF/d®) 2 \ 

G cS = G*A (1 + a ) - - [ 2ZF + 3{dF/m2 ) ■ (3-5) 

The contribution G*A 2 is due to the exchange of a graviton between the two bodies, whereas 
G*A 2 a 2 = G*(dA/d(p) 2 comes from the exchange of a scalar particle between them. Of 
course, when the distance between the bodies becomes larger than the inverse mass of 
the scalar field, its influence becomes negligible and one gets G e R ~ G^. Note that as 
usual, the last expression in Eq. ( p.5| ), in terms of Jordan-frame notation, is much more 
complicated than its Einstein- frame counterpart. In the particular Brans-Dicke representa- 
tion, F = $ and Z = c<j($)/$, it however reduces to the simpler (and well-known) form 
Geff = G^- 1 (2u + 4)/(2lu + 3). 

The experimental bound (|3.2|) shows that the present values of G e s and Gn differ by 
less than 0.02%. However, they can a priori differ significantly in the past. It should 
be noted that the experimental limit on the time variation of the gravitational constant, 
\G e s/G e s\ < 6 x 10~ 12 yr _1 [[40], does not imply any constraint on 2A/A = —F/F. Indeed, 
G e ff can be almost constant even if A (or F) varies significantly. A simple example is provided 
by Barker's theory p3|] , in which A(<p) = cos(p : One gets G e flf = G*(cos 2 <£> + sin 2 (p) = G*, 



which is strictly constant independently of the time variations of A((p(t)). Nevertheless, as 
pointed out in |R9|, under reasonable cosmological assumptions, one can derive G e s ~ G^ 
with ~ 10% accuracy up to redshifts z ~ 1. 



IV. SCALAR-TENSOR COSMOLOGY 

The equations derived in this section generalize those of our previous paper [KBf in several 



ways. First, we use the most general representation ( |2.1| ) of the theory, instead of the simpler 



choice Z = 1 that was made in pB |. Second, we take into account a possible spatial curvature 



of the universe, which will be an interesting possibility in our studies of Sec. VI below. Third, 
we write the equations for an arbitrary pressure of the perfect fluid describing matter in the 
universe. This will not be useful for our reconstruction program of the following sections, as 
matter can be assumed to be simply dustlike for the redshifts z < 5 that we will consider, but 
these general equations may be interesting for further cosmological studies of earlier epochs 
of the universe. Finally, we comment on the Einstein-frame version of these equations, which 
are mathematically simpler, but actually more difficult to use for our purpose. 

A. Background 

We consider a Friedmann-Robertson- Walker (FRW) universe whose background metric 
in the Jordan frame is given by 

ds 2 = -dt 2 + a 2 (t) d£ 2 , (4.1a) 

d f = _^L + r 2 ( d e 2 + sin 2 9 d<j) 2 ) , (4.1b) 

1 — kt v ' 

where k = — 1, 0, or 1 for spatially open, flat, or closed universes respectively. The scalar 
field $ (or tp, in the EF) is also assumed to depend only on time. Since the relation between 
the EF and JF is given by ds 2 = A 2 (if) ds 2 , see Eqs. ( |2.4j ), our universe is still of the FRW 
type in the EF, with dsl = —dt\ + a 2 (£*) dC 2 and 

dt = A((p) dt* , a = A{ip) a* . (4.2) 

In the following, matter will be described by a perfect fluid, and we will write its energy- 
momentum tensor as 

T^ = (p + P)u^u v + pg^ = A' 2 T* v = A' 2 ((p* + p*)w*u* + p*#*J , (4.3) 

where u^ = dx^/\ds\ and u^ = dx^/lds^ are the spacetime components of the four- 
dimensional unit velocity of matter, in JF and EF units respectively. As we are interested 
in a FRW background, the spatial components Ui and u* (i = 1, 2, 3) all vanish. From Q4.3D, 
we deduce the relation between the matter density and pressure in both frames: 

p* = A 4 p, p* = A A p. (4.4) 

The background equations in the JF follow from ( |2.2a| )- (|2.2c|) , and read 

87rG*p + -Z$ 2 - 3HF + U , (4.5a) 

8ttG* (p + p) + Z& + F - H F , (4.5b) 

dF ( ■ 9 k\ dZ $ 2 dU . 

3^r (H + 2H 2 + — -— , 4.5c 

d§ V a 2 J d$ 2 d$ K J 

, (4.5d) 



3F- 


(h 2 + 


a 2 1 


-2F 


■(*- 


a 2 J 


Z 


($ + 3#$) 


P + 3H (p 


+ P) 



where H = d(\na)/dt, and a dot denotes differentiation with respect to the Jordan- frame 
time t. As usual, if p/p = w = const., Eq. (|4.5d|) is trivially integrated as p oc a~ 3 ( 1+w > (and 
in particular p oc a~ 3 for dustlike matter). Equation ( |4.5c| ) is actually a consequence of the 
other three, and we will not need it in the following. 



Since these equations correspond to the most general parametrization ( |2.1| ) of scalar- 
tensor theories, many particular cases are easily recovered. For instance, the case of a 
minimally coupled scalar field |TJ| is obtained for constant values of F and Z (say, F — 1 



and Z = 8irG*), and the particular model considered in ||I] is recovered immediately for 
F = $ and Z = 0. 

The corresponding background equations in the EF are very similar to those in general 
relativity. They follow from Eq. (|2.6a| ), and read 



8vrG>* + bf- + 2V(<p) , (4.6a) 





47rG„(p„ + 3p„) + 2Nf -2V(<p), (4.6b) 



where H* = d(lna*)/dt* is the Einstein-frame Hubble parameter. It is obvious from (|4.6b|) 
that a vanishing potential V(ip) implies d 2 a*/dtl < 0, so that the universe is decelerating 
in the Einstein frame. However, because of the relation a = A(ip) a*, see Eq. ( f4.2|) , the 
observed (Jordan-frame) expansion rate a may be positive even in this case, and we will see 
concrete examples in Sec. VIA below. This is an important point to remember: Although 
we are looking for cosmological FRW backgrounds whose expansion is accelerating, the sign 
of d 2 a*/dtl is a priori not fixed. 

The scalar- field equation of motion in the EF follows from Eq. ( |2.6b| ), and reads 

dv + 3H/ £ + m^ = _ 

dt% dt* dip 

It is also similar to the usual Klein-Gordon equation, with the notable difference of a source 
term on the right-hand side, with the coupling strength at((p) defined in Eq. fl2.7|) above. 

It is tempting to tackle our problem in the EF as the equations are simpler and we can 
rely on experience gained in general relativity. However, a crucial difficulty that we encounter 
is that all physical quantities which appear in the EF background equations are not those 
that come from observations. Moreover, the behavior of matter in the EF is complicated by 
the relations ( [4.41 ): Instead of the simple power law p oc a -3 for dustlike matter in the JF, 
one gets p* = A 4 p oc Aa~ 3 in the EF, where A(cp(aA) can have a priori any shape. To avoid 
these problems, we will thus work in the JF, and show that the "reconstruction" program 
can equally well be implemented, like in general relativity, although it is mathematically very 
different. We will nevertheless check at the end the consistency of the solutions obtained by 
translating them in terms of EF quantities. 

B. Perturbations 

We now consider the perturbations in the longitudinal gauge. For this problem, we will 
restrict our discussion to the case of a spatially flat FRW universe (k — 0), and write the 
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JF and EF metrics as 



ds 2 
ds 2 = 



-{l + 2<j))dt 2 + a 2 {l-2?lj)dx 2 , (4.8a) 

1 + 2(j)*)dt 2 + a 2 (l - 2^)dx 2 . (4.8b) 

In the EF, the perturbation equations deriving from Eq. ( |2.6aj ) are strictly the same as in 
general relativity plus a minimally coupled scalar field. One thus finds notably 0* =ip*. On 
the other hand, the equations for scalar-field and matter perturbations are modified by the 
matter-scalar coupling, proportional to oc((p) in Eqs. ( [2.6b|) and (|2.6c|) . 

For our purpose, it will be more useful to write the perturbation equations in the (phys- 
ical) JF. Let us define the gauge invariant quantity! 

6p 



Sn 



P + P 



+ 3Hv 



(4.9) 



where v is the matter peculiar velocity potential (such that 5u^ = —d^v is the perturbation 
of the four-dimensional unit velocity u^). We now work in Fourier space, and assume a 
spatial dependence exp(zk ■ x), with k = |k|. The conservation equations of matter Q2.2c|) 
give 



k 2 d(tp + Hv) 

-—v + 3- 



dt 



P 



v + - (2Hv - 8 ri 
P 



(4.10a) 
(4.10b) 



On the other hand, the Einstein equations ( |2.2a| ) give 

i/j = + 5F/F , 
2F(iP + H(f>) + Fcf) = MG,(p + p)v + Z$6$ + 6F-H6F 



(4.11a) 
(4.11b) 



-3F0 - [ 2— F - Z& + 3HF 



7 j = 8ttG # (p + p)5 m + f ^ - 6H 2 -3^\SF + SU 



a' 



F 2 



+Z$<J$ + 3HZ$ 5$ + -5Z $ 2 + 3- SF 



(4.11c) 



Note that ^ ip in the JF, in contrast to the corresponding problem in general relativity 
or in the EF. Equation ( 4.1 la|) is actually an obvious consequence of the relation between 
g^ u and g^ u , Eq. fl2.4a| ), and of the fact that 0* = ■0*. Finally, Eq. ( |2.2b| ) yields the equation 
for the dilaton fluctuations 5$: 

sq+Uh + ^^s* 



+ 



K-m + 2H 2 )±( l - d -f\ + ±( l - d ^ 



k- 



d$ \Z d$J 



rf$ \Z d$ 

1 dF 
Zd$ 



+ 



d 2 lnZ $' 



d<& 2 

+ (3V> + 0)6 - 



5$ 



dU 
Zd$ 



(4.12) 



3 Note that our definition differs from the quantity e m introduced in 



: e r 



(l+p/p)5„ 
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In the particular representation Z = 1 used in Ref . |36| , this equation reduces to the simpler 
form 



<5$ + 3#<5$ 



5$ 



k 2 • TT2 ,d 2 F d 2 U 

a 2 gM> 2 d$ 2 

k 2 

— (0 - 2^) - 3(V> + 4#V + #</>) 



_ + (3 ,) + ,)4-2^-. (4.13) 



V. THE RECONSTRUCTION PROBLEM 



The reconstruction of the potential U(Q) was shown in |14]] to be possible in the frame- 
work of general relativity plus a minimally coupled scalar field, the A-field or quintessence, 
provided the Hubble diagram (and thus also H(z)) can be extracted from the observations. 
An essential difference arises when one deals with scalar-tensor theories: We have to recon- 
struct two unknown functions instead of one, hence we need to extract two quantities (as 
functions of the redshift z = a /a — 1) from the observations. Actually, in the minimally 
coupled case, the knowledge of the luminosity distance D L and of the clustering of matter 
5 m , both in function of z, provides two independent ways to reconstruct the scalar field 
potential [jyJj.H In our case, both quantities are necessary and the reconstruction itself is 
significantly more complicated. 

The present section generalizes our previous results of Ref. |36| not only by considering 



the most general parametrization ( |2.1| ) of scalar-tensor theories and by taking into account 
the possible spatial curvature of the universe, but also by discussing particular cases that 
were excluded in this reference. From now on, we will restrict our discussion to the case of a 
pressureless perfect fluid {p = = p*), because all matter in the universe will be assumed to 
be simply dustlike, of course besides that part needed to account for the present accelerated 
expansion (i.e., the scalar field in the present framework). 



A. Background 

The first step of the reconstruction program is the same as in general relativity, since it is 
purely kinematical and does not depend on the field content of the theory: If the luminosity 
distance D L is experimentally determined as a function of the redshift z, one can deduce 
the quantity H (z) from the relation 



4 More precisely, to reconstruct the potential U(&) without any ambiguity in the minimally coupled 
case, one needs to know both Dl(z) and the present energy density of dustlike matter £l m ,o-, or 
both 5 m (z) and the present value of the Hubble constant Hq. In ° ur general scalar-tensor case, we 
need to know the two functions Dl(z) and 5 m (z), but no independent measurement of Q m fi or Hq 
is necessary. 
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Hz) 




i + n 



'H D L (z) 



K,0 



-1/2 



(5.1) 



where the prime denotes the derivative with respect to z. The large square brackets contain 
a corrective factor involving the present energy contribution f2 Kj0 = — K /( a o-^o) of the spatial 



curvature of the universe. It was not written explicitly in Refs. [||,|14]], which focused their 
discussions on the flat-space case (fi Kj o = 0), but it is a straightforward consequence of 
Eqs. (23)-(25) of Ref. ||. Since present experimental data suggest that \Q K> o\ is small, the 
flat-space expression for 1/H(z) = [D L (z)/(l + z)]' is a priori a good approximation anyway. 
Note that even if one uses the exact equation fl5.1|), it reduces to the flat-space expression 
for z = (because -Dl(O) = 0), and therefore ifo is always known without any ambiguity. 
To determine H(z) precisely at higher z, one then needs to know both Dl(z) and Q K fi- 

By eliminating Z$ 2 from the background equations (|4.5a|) and (|4.5b|) , we then obtain 
the equation 



F + 5HF + 2(H + 3H 2 + -^jF = 8nG*p + 2U , 
which, when rewritten in terms of the redshift z, gives the fundamental equation 



(5.2) 



F" + 



(In if)' 



4 



1 + z 



F' + 



6 



1 + z) 



T^^'- 4 (f)^ 



217 



:i + ^) 2 ff 2 



+ 3(1 + z) 



Ho 
H 



fb^m.O 



(5.3) 



As before, an index means the present value of the corresponding quantity, and we use again 
the notation /' = df/dz. In this equation, f2 m0 = 87tG*Po/(3F q Hq) stands for the present 
energy density of dustlike matter relative to the critical density e crit = 3ifQ/87rGjv,o- To 
simplify, this critical density is defined in terms of the present value of Newton's gravitational 
constant (|3.4|) , Gjv,o — G*/F , instead of the effective gravitational constant ( |3.5D actually 
measured in Cavendish-type experiments. Indeed, solar-system experiments tell us that 
their present values differ by less than 0.02%, as discussed in Sec. III. [Note in passing that 
by changing the value of G>, one can always set F = 1 without loss of generality] 

In conclusion, we are left with a non-homogeneous second order differential equation 
for the function F(z), a situation very different from that prevailing in general relativity. 
However, the right-hand side also depends on the unknown potential U(z), so that this 
equation does not suffice to fully reconstruct the microscopic Lagrangian of the theory. 
As we will show in Sec. VI below, it can nevertheless be used for a systematic study of 
several scalar-tensor models, provided one of the two unknown functions is given (or a 
functional dependence between them is assumed). This can be useful as we do not expect a 
simultaneous release of data yielding H(z) and S m (z). We will see that such a study already 
yields powerful constraints on the family of theories which are viable. 

On the other hand, if S m (z) is also experimentally determined, and if we assume a 
spatially flat FRW universe (£l K = 0), we will see in the next subsection (V.B) that the 
value of fl m fi as well as the function F(z) can be obtained independently of U(z). Equation 
then gives U(z) in an algebraic way from our knowledge of H(z), F(z) and fl m fi. 
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Let us now assume that both F(z) and U(z) are known, either because one of them was 
given from theoretical naturalness assumptions, or because S m (z) has been experimentally 
determined with sufficient accuracy. We will also assume that both f2 m0 and f2 Kj0 are 
known. It is then straightforward to reconstruct the various functions of $ entering the 
microscopic Lagrangian ( |2.1| ). In the Brans-Dicke representation, one has F = <&, therefore 
the knowledge of F(z) and U(z) suffices to reconstruct the potential U(<&) in a parametric 
way. However, to fully determine the theory, one also needs to know u>(<&) = $Z($), or 
equivalently an equation giving the z-dependence of Z. On the other hand, in the simpler 
representation Z = 1 and F(<&) unknown, we need an equation giving the z-dependence 
of $ to reconstruct F(3>) and U(<&) parametrically These two cases, as well as any other 
possible parametrization of the theory, are solved thanks to Eq. ( |4.5b|) above, which reads 
in function of the redshift 



Z $' 



or equivalently 



1 



Z $ 



/2 



3F' 



+ 



-F" 

-3(l + z 

3F 



(ln#)' + 



l + z 



F' + 2 



(InH)' 

l + z 



IT) n *-° 



^) FaQ 



l + z [l + zf 



H 



3F 



i4 m,0 5 



(5.4) 



H 



^K,0 ~~ 



U 



'l + z) 2 H 2 



-3(l + z) 



H, 



H ) 



F n Q 



0"m,0 • 



(5.5) 



In the Z = 1 representation, §(z) — $ is thus obtained by a simple integration. In the 
Brans-Dicke representation, on the other hand, u)(z) is given by an algebraic equation in 
terms of H(z), F(z) = $(z), and their derivatives. 

It is rather obvious but anyway important to note that if the microscopic Lagrangian 
( ^■1| ) can be reconstructed in the JF, it can also be directly obtained in the EF, Eq. ([2.5|). 
This allows us to check the mathematical consistency of the theory, and notably if the 
helicity-0 degree of freedom tp always carries positive energy. One can also verify that the 
function A(ip) defining the coupling of matter to the scalar field is well defined, and notably 
single valued. Finally, the second derivative of the potential V(ip) also gives us the sign of 
the square of the scalar mass, and negative values would strongly indicate an instability of 
the model. These important features cannot easily be checked in the JF, because the sign 
of Z($) in Eq. fl2.1p is not directly related to the positivity of the scalar-field energy (see 
below), and also because the second derivative of U(<&) does not give the precise value of 
its squared mass. [As shown by Eq. ( p.4d|) , the helicity-0 degree of freedom tp may have a 
mass, d 2 V (ip) / dtp 2 ^ 0, even if £/($) is strictly constant, provided F{$F) varies.] 

Let us thus assume that H(z), Q m> o and Q Kt0 are known, and that F(z) and U(z) were 
reconstructed as above. Equation ( 2.4c|) then gives A(z) = -F~ 1//2 (z), i.e., the Einstein- 
frame coupling factor A as a function of the Jordan-frame redshift z (which is the redshift 
we observe). Combining now Eq. ( f2.4b| ) with ( |5.4| ), we get 





+ 



Z$ 



f'2 



IF 



(5.6a) 
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3 (F'\ 2 F" 
~ 4 \F) 2F 


\ 1 (\nH)'+ 1 

L2 V ; l + z_ 


F' (In if)' 

F + l + z 


~\h) K '°~2~(' I \W) y fim - ' 


or also, not eliminating the potential U 


(d<pY_3(F>y 3F> 3 O ^o\ 2 n 
\dz) ~ A\f) (1 + z)F ' (1 + *)' 6 {h) "*° 


u 


3(l + .)(f) 2 f^,o. 


(1 + z) 2 FH 2 



(5.6b) 



(5.7) 

The EF scalar ip is thus also known as a function of the Jordan-frame redshift z (up to an 
additive constant <p which can be chosen to vanish without loss of generality), and one can 
reconstruct A((p) in a parametric way. Similarly, the EF potential V(<p), Eq. fl2.4d| ), can be 
reconstructed from our knowledge of F(z), U(z) and <f(z). 

Since ip describes the actual helicity-0 degree of freedom of the theory, this field must 
carry only positive energy excitations, and {dip/dz) 2 must be positive. On the other hand, 
the tensor and scalar degrees of freedom are mixed in the JF, and the positivity of energy 
does not imply that ZQ' 2 should always be positive. Actually, Eq. (|5.6a| ) shows that it can 
become negative when |(lnF)' 2 happens to be larger than ip' 2 , which can occur in perfectly 
regular situations. [We will see an explicit example in Sec. VIA below.] This underlines that 
the parametrization Z = 1 can sometimes be singular: The derivatives of $ may become 
purely imaginary although the scalar degree of freedom ip is well defined. On the other 
hand, the Brans-Dicke representation is well behaved (<J> /2 remains always positive), and the 
positivity of energy simply implies the well-known inequality u>(3>) > — |. Actually, the 
particular value u = — | is also singular, as it corresponds to an infinite coupling strength 
a = (2u + 3) -1//2 between matter and the helicity-0 degree of freedom ip. The domain for 
which the Z = 1 parametrization is pathological although the theory remains consistent 
simply corresponds to — | < tu(&) < 0, or \a\ > l/\/3. 



B. Perturbations 

Although the perturbations will not be used in Sec. VI below, we emphasize that the 
phenomenological reconstruction of the full microscopic Lagrangian can be implemented 
without any ambiguity if fluctuations are taken into account. For completeness, we review 
now this part of our program. We assume that both H(z) and the matter density perturba- 
tion S m (z) are experimentally determined with enough accuracy, and as in Sec. IV.B above, 
we focus our discussion on the case of a spatially flat FRW universe (0 K = 0). We also 
assume that matter is dustlike (p = 0), and the perturbation equations of Sec. IV.B are thus 
simplified. In particular, Eq. (|4.10b] ) reduces to the mere identity <ft = v. 

We consider comoving wavelengths A = a/k much shorter (for recent times) than the 
Hubble radius H , and also shorter than the inverse mass of the scalar field: 

k 2 /a 2 ^>max[H 2 ,A- 2 \d 2 V/dip 2 \) . (5.8) 
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Two different reasonings can now be used to reach the same conclusions. The first one, 



explained in Ref . |36| , consists in taking the formal limit k — > oo in the various perturbation 
equations. Then, the leading terms are either those containing d~ m or those multiplied by the 
large factor k 2 /a 2 . One also needs to consider only the growing adiabatic mode of Eq. ( (4.12|) , 
for which \5&\ < fc 2 a" 2 |<5<l>|. 

The other reasoning needs a simpler (but a priori stronger) hypothesis. One assumes 
that the logarithmic time derivative of any quantity, say /, is at most of order H : |/| < \Hf\. 
Physically, this means that the expansion of the universe is driving the time evolution of 
every physical quantity. Then the hypothesis k 2 /a 2 ^> H 2 suffices to derive straightforwardly 
all the following approximations. 

Note that both reasonings correspond in fact to the same physical situation of a weakly- 
coupled light scalar field. In the case of a strongly-coupled but very massive scalar (see the 
second paragraph of Sec. Ill), the equations cannot be approximated as shown below, and 
the time evolution of density fluctuations does not follow the same law. For instance, in the 



particular model considered in Ref. |[41|| , one always finds a strong clustering of the scalar 
field at small scales. Indeed, this model corresponds to the choice F = $ and Z = in 
action (pT|) , and Eq. ( gTg) can then be rewritten as (d 2 U/d® 2 )5$> = (fc 2 /a 2 )(0-2^)-3(^ + 
AHtp + H<p) — 2<f)(dU / d&) . Therefore, even if the scalar field is very massive (d 2 U/d& 2 large), 
one finds that it is anyway strongly clustered for comoving wavelengths a/k shorter than 
the inverse mass, i.e., in the formal limit k —>■ oo. Although this is a priori not forbidden by 
observations of gravitational clustering, since the inverse mass must be much smaller than 
the astronomical unit in this model, this is anyway an indication of its probable instability. 
We will not consider such heavy scalar fields any longer in this paper, and we now come back 
to the class of weakly-coupled light-scalar models, which are the most natural alternatives 
to general relativity. 



Setting B = ip + Hv and making use of ( |4.10b|) , one can write (|4.10aj) as 



k 2 
5 m + 2H5 m + — = W + 6HB « , (5.9) 



a? 



where the right-hand side is negligible with respect to each separate term of the left-hand 
side because of the above hypotheses. Note that ( |5.9j ) just reproduces the standard evolution 
equation for matter perturbations. Using (|4.13|) , we also arrive at 

M. rj. n,-, dF / d ^ j. FdF/d<S> 

M,a ^- 2 * ) -Z- a -* ZF + UdF/M)' ' (5 ' 10) 

where the second equality is a consequence of Eq. (|4.11a| ) . In the case of GR plus a minimally 
coupled scalar field, one finds that 5$ oc k~ 2 <p in the limit k — > oo, so that the scalar field is 



not gravitationally clustered at small scales [14 1. This is in agreement with the observational 



fact that the dark matter described by the A-term should remain unclustered up to comoving 
scales R ~ 10 h^^+z)' 1 Mpc (where we recall that h' 1 = 100 Hq 1 km s" 1 Mpc -1 ). On the 
other hand, in our scalar-tensor framework, Eq. ( |5.10|) shows that the scalar field is clustered 



at arbitrarily small scales, but only weakly because the derivative \dF/d$\ is experimentally 
known to be small [see the solar-system constraint ( |3.2|) , and the limit a 2 ^$ 0.1 justified 
5B| for redshifts z ^ 1] . The class of models we are considering, involving a light scalar 



in 
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field weakly coupled to matter, is thus also in agreement with observations of gravitational 
clustering. 

Finally, still under the above hypotheses, Eq. ( |4.11c| ) implies 

k 2 k 2 dF 

-2—F ( j ) n8TrG*p6 rn + ——6®. (5.11) 

a z a 2 d<P 



Remembering the definition ( |3.5| ) for G e g, and using ( |5.10| ) above, Eq. ( |5.11| ) can be recast 
in a form which exhibits its physical content: 

^<j> « -AnG cS p5 m . (5.12) 

a 2 

Poisson's equation is thus simply modified by the substitution of Newton's constant G by 
Geff, the effective gravitational constant between two close test masses! This conclusion was 



also reached in ||3q] , but only for Brans-Dicke theory with a constant parameter u, while 
we have derived it for an arbitrary (light) scalar-tensor theory. As discussed in Sec. Ill 
above, expression ( p~5| ) is valid only if the distance between the test masses is negligible 
with respect to the inverse scalar mass. The physical reason why this expression appears 
in Poisson's equation ( |5.12| ) is just that we are working in the short wavelength limit 



The frequency of the waves we are considering is so large that the scalar field behaves as if 
it were massless. 

Combining fl5.9| ) with (|5.12 ), we now arrive at our final evolution equation for S m : 



5 m + 2H5 m - 4nG cS p 5 m « . (5.13) 

In terms of the redshift z, this reads 

H 2 C + (&L - J^j 8> m « \{1 + z)Hl^0± Qmfi 6m . (5.14) 

Provided we can extract from observation both physical quantities H(z) and S m (z) with 
sufficient accuracy, the explicit reconstruction of the microscopic Lagrangian is obtained in 



the following way. Starting from (|5.14f ) and using the fact that today GW,o and Gn,o differ 
by less than 0.02%, Eq. ( |5.14| ) evaluated at present gives us the cosmological parameter Q m ,o 
with the same accuracy. Then, returning to Eq. (|5.14j) for arbitrary z, we get G e s(z) = p(z), 
where p(z) is a known function of the observables H(z), 5 m (z), and their derivatives. Using 
now Eq. (|5.4fj and expression ( |3.5|) for G cff , we get a nonlinear second order differential 
equation for F(z), which can be solved for given F and Fq [one can always set F = 1 
without loss of generality, while Fq is constrained by Eq. ( |3.2j )1. After we have found F(z), 



we can plug it into ( [5 .3D to determine U(z) in an algebraic way. The final step is explained 



in the previous subsection, above Eq. ( |5.4j ), for the various possible parametrizations of 
action (|2.1| ): In the Z — 1 parametrization, $(z) — $o is obtained by a simple integration of 
Eq. ( |5.4|) , while in the Brans-Dicke parametrization (F(<&) = $), u>(z) is given algebraically 
by the same Eq. (|5.4j) . This enables us to reconstruct F(§) (or uj(§)) and U($>) as functions 
of $ — $o for that range corresponding to the data. 

Actually, for sufficiently low redshifts z ^ 1, Eq. ( |5.14| ) can be simplified without losing 



too much accuracy. Indeed, as shown in Ref. []55| , the square of the matter-scalar coupling 
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strength a, Eq. fl2.71) , is at most of order 10% for such redshifts. Moreover, under natural 
assumptions, much smaller values of a 2 are generically predicted in scalar-tensor theories 
|25| , |26|| . Therefore, G e fj and GV differ by less than ~ 10% for redshifts z < 1, and Eq. ( [5. 14| ) 



can be used to obtain G e g/GN,o ~ Gn/Gn,o — F /F with the same accuracy. The interest 
of this simplification is that F(z) is now given by an algebraic equation. In the Brans-Dicke 
representation, all the steps of the reconstruction program are thus algebraic, Eq. ( |5.3| ) 
giving U(z), and Eq. ( |5.4| ) giving uj{z). The only non-algebraic step is the final parametric 
reconstruction of £/($) and a;($). 

Let us end this section by a few comments on the observational accuracy which will 
be needed for this reconstruction program to be implemented. First, Eq. ( |5.14j ) allows to 
reconstruct F(z) only if S' m and 5'^ are both determined with enough accuracy. Moreover, the 
second derivative of this reconstructed F(z) is needed in Eq. ( |5T3| ) to obtain U(z). Therefore, 
the actual reconstruction of the potential depends a priori on the fourth derivative of S m (z), 
so that extremely clean data seem to be necessary. However, the situation is better than this 
naive derivative counting suggests. Indeed, the above estimates for a 2 show that F(z) does 
not vary much on the redshift interval < z ^ 1. Therefore, the first two terms of Eq. ( p.3| ), 
involving F' and F", are expected to be negligible with respect to the third one involving 
F. A noisy experimental determination of 6'^(z) and 8'™(z) is thus not a serious difficulty 
for our reconstruction program. On the other hand, clean enough data are still needed to 
determine F(z) from Eq. ( |5.14| ), using 5 m (z) and its first two derivatives. Before such clean 
data are available, it will be sufficient to verify that Eq. ( |5.14j ) is consistent with a slowly 
varying F(z). In the next section, we will show that interesting theoretical constraints can 
anyway be obtained without knowing at all the density fluctuation 5 m (z), but using only 
the luminosity distance D L (z) and consistency arguments within particular subclasses of 
scalar-tensor models. 



VI. CONSTRAINTS FROM AN ACCELERATING UNIVERSE 



In Ref. 0], a fit of presently known supernovae events has been performed to obtain 
the luminosity distance Dl(z) up to redshifts z ~ 1, of course still with large uncertainties. 
Although this is not yet sufficient to constrain seriously scalar-tensor models, we can expect 
clean data on Dl(z) in the near future from additional supernovae events, and anyway 
earlier than for the density perturbations 5 m (z). The SNAP satellite will in particular 
observe thousands supernovae events up to z ~ 1.7. In this section, we will concentrate on 
the theoretical constraints that can be extracted from the knowledge of Dl(z) alone, and 
therefore of H (z) using Eq. ([5.11). We will thus only use the results of subsection V.A above. 
Since the knowledge of this function does not suffice to fully reconstruct the microscopic 
Lagrangian (|2.1| ), we will need additional assumptions on one of the functions it involves, 
either F (or Z, depending on the parametrization) or the potential U. One may also assume 
a functional relation between F and U (for instance U oc F M as in Ref. f45fl). 



To emphasize as clearly as possible what kind of constraints can be imposed on scalar- 
tensor theories, we shall consider the worst situation for them. Let us assume that the 
observed function H(z) will be exactly given by Eq. ( |4.5a| ) for k — 0, F — 3> = 1, and 
U = A = 3H 2 n Afi : 



(H/H y = fi A ,o + On,o(l + zf . (6.1) 

Of course, such an observation would a priori call for the following standard interpretation: 
Gravity is correctly described by general relativity, and we live in a flat universe filled with 
dustlike matter and a cosmological constant, with corresponding present energy densities 
(relative to the critical density) Q m ,o and fl\,o- However, for our purpose, Eq. ( |6.1| ) should 
just be considered as kinematical. It tells us how the universe expands with redshift z, but 
we are free to assume that the dynamics of the expansion is governed by a scalar-tensor 
theory. Therefore, £l m fi and Oa,o are here mere parameters, whose names refer to their 
physical significance in the framework of GR. Of course, one should not forget that they do 
not have the same interpretation within scalar-tensor theories. 

For our numerical applications, we will further take the present estimates based on 
combined CMB fluctuations and supernovae observations (they will be determined more 
accurately by future experiments): 

Q Afi ~ 0.7 , On.o « 0.3 . (6.2) 



For these numerical values, fl6.1|) is consistent with the presently available luminosity distance 
Dl(z) up to z ~ 1. Actually the best-fit universe, if we assume flatness, gives f2\.,o = 0.72 
and fl m ,o = 0.28. We have chosen to work directly with the exact form (|6.1|), instead of 
the Dl(z) extracted from observation, in order to clarify the physical content of our results. 
Indeed, the present observational estimates for Dl(z) are still too imprecise to constrain 
strongly the class of scalar-tensor theories we are considering. Moreover, some of our results 
below depend crucially on the fact that H(z) keeps the form fl6.1| ) up to redshifts z ~ 2, 
which have not yet been reached experimentally. To relate our results to those obtained in 
T6||i6|| using fitting functions or an expansion in powers of z, one just needs to use Eq. fl5.1|): 



Our exact expression (|6.1|) for H(z) corresponds to some exact expression for D L (z 

To summarize, we are assuming in this section that future observations of the luminosity 
distance Dl(z) will provide a H(z) of the form ( |6.1|) with the numerical values (|6.2|) . This 
implies notably that our Universe is presently accelerating. On the other hand, we are not 
assuming that the correct theory of gravity is necessarily GR plus a cosmological constant. 
The main question that we will address is therefore the following: Would such an "observed" 
H(z) necessarily rule out the existence of a scalar partner to the graviton? If not, would 
it be possible to reproduce ( |6.1|) within a more natural scalar-tensor theory, in which f^o 
could be explained by a "generalized quintessence" mechanism? 

We will first analyze in subsection A the simplest subclass of scalar-tensor theories that 
we can consider, namely when U = in action ( |2.1| ). Since this is a priori the subclass 
which differs the most from GR plus a cosmological constant, this study will be rather 
detailed, and it will allow us to underline the mathematical and physical meaning of the 
constraints that are obtained. Subsection B will be again devoted to the case of a massless 
scalar field, but combined with a cosmological constant. As its conclusions basically confirm 
those of subsection A, we will present them more concisely. Finally, we will briefly discuss 
in subsection C the cases where one imposes particular forms for the coupling function F in 
action ( j2.1| ), and one reconstructs the potential U from the background equations ( |5.3| )-(3~E). 
The case of a given functional dependence between F and U will also be addressed. 
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A. Case of a vanishing scalar-field potential 

Since a cosmological constant can be interpreted as a particular case of scalar- field poten- 
tial, it is instructive to analyze whether an observed expansion like (|6.1|) could be reproduced 



in a theory without any potential, and we now study Eqs. (|5.3|) - (|5.6|) for £/($) = = V{ip). 
This case can be analyzed using the second order differential equation ( |5.3| ) for F, which 
simplifies significantly if one introduces a function / such that 

F(z)/F = (l + z) 2 f(l + z). (6.3) 

[As mentioned in Sec. V.A above, one can also set F = 1 without loss of generality] Then, 
using the assumed "experimental" expression ( |6.1| ) for H(z), and writing ( |5.3| ) in terms of 
x = 1 + z, we get 

(O A)0 + n m ,o x 3 )xf"(x) + hl mfi x 3 f\x) - 4Q Kfi xf(x) = 3Q mfi . (6.4) 

To avoid any confusion, let us recall that Qa,o (and the two occurrences of fl m ,o i n the 
left-hand side) comes from the "observed" cosmological function ( |6.1| ), notwithstanding the 
fact that there is no cosmological constant in the model we are considering. The value Q m ,o 
appearing in the right-hand side stands for the present relative energy density of dustlike 
matter. We assume that it takes the same numerical value (|6.2|) as in the "observed" H(z) 
( p.lp . Equation ( |6.4| ) tells us how we should choose f(x) to mimic exactly this H(z) in the 
present potential-free theory. In other words, ^a.o an d ^ m ,o are t wo numbers assumed to be 
given by experiment, and we wish to fit f(x) and Q Kt0 to satisfy Eq. Q6.4J ). 

To integrate this second-order differential equation, we need two initial conditions for / 
and its derivative. The first one is an obvious consequence of Eq. (|6.3|) taken at z — 0, and 
we simply get /(l) = 1. The second one should be such that the solar-system bound ( |3.2| ) is 
satisfied. For instance, if $[, does not vanish, it is sufficient to impose Fq = 0, i.e., f'{l) = —2 
using Eq. ( |6.3p . This corresponds to a scalar-tensor theory which has been attracted towards 
an extremum of F during the cosmological expansion of the universe (cf. p5|j26[1 ), so that 



it is presently strictly indistinguishable from general relativity in solar-system experiments. 
[The full allowed domain for f'(l) will be explored below in a numerical way] 

1. Spatially flat universe 

We consider first our potential- free model in a spatially flat FRW universe (Q Ky o = 0). 
Then Eq. (|6.4j) becomes a first-order differential equation for /', and its integration yields 



f'(x) 



VTTCx-z 



yi + (x*-i 




2J1 + C 



(6.5) 



where we have set ( = Q mt0 /Q^o, and where the final constant inside the square brackets 
has been chosen to impose /'(l) = —2 (i.e., Fq = 0). The function f(x) = l + f*f'(y)dy can 
be explicitly written in terms of generalized hypergeometric functions, but its complicated 
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expression will not be useful for our purpose. Let us just quote the first order of its expansion 
in powers of fi m ,o/^A,o : 



/(*) 



2x + - (15 - 16x + x 4 + 12z In x) -^ + O 
4 V '^a.o 



^m,0 
< ^A,0 , 



(6.6) 



In conclusion, Eq. (|6.4j ) could be integrated analytically, in the particular case of a spatially 
flat universe. This means that at least in the vicinity of z — 0, there a priori exists a 
potential-free scalar-tensor theory which exactly mimics general relativity plus a cosmolog- 
ical constant. 

However, the theory is mathematically consistent only if F(z) remains strictly positive. 
[If F vanishes, then the coupling function A((p), Eq. ( |2.4cj ), between matter and the helicity- 
degree of freedom ip diverges, and if F becomes negative, the graviton carries negative 



energy] Let us thus compute the value z, 



mnx for which F(z max ), 



or /(l + z max ), vanishes 



for the first time. Because of the complexity of the solution /, we did not find a close 
analytical expression for z ma , x , but its expansion in powers of fi m ,o/^A,o can be obtained 
straightforwardly: 
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(6.7) 



Numerically, for the values (|6.2j ) of Q\ t0 and fi m0 , we find z max ~ 0.66. In conclusion, this 
scalar-tensor model is able to mimic general relativity plus a cosmo logical constant, but only 
on the small interval z < 0.66. If future observations of type la supernovae give a behavior 
of H(z) of the form ( |6.1| ) on a larger interval, say up to z ~ 1, then the present scalar-tensor 
theory will be ruled out. This example of a vanishing potential illustrates a conclusion that 
we will reobtain below for more general theories: The determination of the form of H(z) 
over some (even rather small) redshift interval is in fact more constraining than the precise 
value of the parameters Q m ,o, ^a,o themselves. Indeed, Eq. (|6.7f ) clearly shows that z max 
cannot exceed 1 even in the presumably unrealistic case of Q m fi ~ ^a,o- [A calculation using 
the exact expression for f(x) shows that z max would exceed 1 only for f2 m0 /^A,o > 1.59.] 
Note that all the results obtained are independent of the parameter H . 



2. Spatially curved universe 

One could try to increase z max by considering a spatially curved FRW universe. We did 
not solve Eq. (|6.4| ) in the most general case, but since we wish to compute the corrections 
to Eq. ( |6.7D due to a small value of \fi K ^/flA,o\, h is sufficient to work at zeroth order in 
^m,o/^A,o- Let us thus set Q mt o = in Eq. ( |6.4|) , which reduces to 



n A ,of"(x)-4:Q K ,of(x) 







(6.8) 



Its solution is obviously a sine if fi Kj0 < (i.e., k = +1, closed universe), or a hyperbolic 
sine for Q K} q > (i.e., k = —1, open universe). Taking into account the initial conditions 
f(l) = 1 and /' (1) = -2, we thus get 
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f(l + z)=co8(2Sz)-\sm{2Zz) for £ 2 = -^° > , (6.9a) 

? "A,0 

/(l + z) = cosh(2£s) - 7 sinh(2£z) for £ 2 = +^ > . (6.9b) 

The first zero of /(l + z) is then reached either at z max = ^ arctan£ or at -k arctanh£. In 
both cases, the expansion in powers of £ gives z max « ~+~ ^ k ,o/^a,o- Working perturbatively, 
one can also compute the correction to this expression due to the nonzero value of fi m> o> and 
one finds that z max is given by Eq. ( |6.7|) above plus the following correction: 



Sz = ± 
6 



12 



l + 4561n(3/2)-163 ^o + / ^,o 

16 fi A , V^AO, 



2 



^ +0 (k)- (610) 



In conclusion, z m&K can be slightly enlarged if we consider our potential-free scalar-tensor 
theory in an open FRW universe (fi Ki0 > 0). Numerically, for the values ( |6.2|) of fi\ t0 
and fl m ,o, we find <5z max ~ 0.26f2 Ki0 /fiA,o- Since the latest experimental data on CMB 
temperature fluctuations already constrain \Q K) o\ to be small (see the latest Boomerang and 
Maxima data), and actually an open universe is unlikely while a marginally closed universe is 
still acceptable, we thus recover the same qualitative conclusion as in the spatially flat case: 
It is possible to mimic general relativity plus a cosmological constant within a potential-free 
scalar-tensor theory only on a small redshift interval z < 0.8. 

3. Numerical integrations 

The above conclusions have been confirmed by numerical integrations of Eqs. fl5.3|) -( |OD , 
still assuming a Hubble diagram consistent with (|6.1|) . Instead of considering only theories 
which are presently indistinguishable from general relativity (Fq = 0), we imposed arbitrary 
initial conditions for F', and computed the corresponding value of the present scalar-matter 
coupling strength a , Eq. (|2.7|). In the case of a spatially flat FRW universe, we recovered 
that the solar-system bound ( |3.2j ) imposes the limit z max ~ 0.68, consistently with the above 
analytical estimate ( |6.7| ). In other words, the constraint ( |3.2[ ) is so tight that even taking 
the largest allowed value for |ao| does not change significantly z max . Figure 1 displays the 
reconstructed F(z) for this maximal |«o|, and one can note that its slope at z = is visually 
indistinguishable from the horizontal. This figure also plots the Einstein-frame scalar ip, 
Eq. (|2.4b| ), which is the actual helicity-0 degree of freedom of the theory. Notice that it 



diverges at z mSLX , so that the theory loses its consistency beyond this value of the redshift. 

Curiously, we found that even if no experimental constraint like (|3.2|) is imposed on |cto| 
(i.e., even if we forget that solar-system experiments confirm very well general relativity), 
then the mathematical consistency of the theory anyway imposes z < 3.5. In fact, Eq. ( |5.3| ) 
alone can be solved for arbitrary large values of z, i.e., there exist initial values of Fq such 
that F(z) remains positive for any z. However, the values of Fq needed to integrate Eq. ( |5.3| ) 
beyond z = 3.5 correspond to negative values of a% = (2u + 3) _1 (where u denotes the 
present value of the Brans-Dicke parameter). In other words, the expression of (dip/dz) 2 
given by Eq. (|5.6j ) would become negative around z = 0, and the helicity-0 degree of freedom 
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would thus need to carry negative energy at least on a finite interval of z, if one wished to 
integrate Eqs. Q5.3|) -( fT6l) beyond z = 3.5. 

Figure 2 displays the maximum redshift z max consistent with the positivity of energy 
of both the graviton and the scalar field, but for any value of the present matter-scalar 
coupling strength \ocq\. As underlined above, one finds that z max can never be larger than 
3.5. This figure also indicates the present solar system bound on |ao|, corresponding to 
Zmax ~ 0.68 as in Fig. 1. The limiting case of a vanishing |ao|, i-e., of a scalar-tensor theory 
which is presently strictly indistinguishable from GR in the solar system, corresponds to 
z max ~ 0.66, as was derived analytically in Eq. (|6.7|). Figure 2 also indicates the range of 



values for |ct | that are generically obtained in Refs. |25[ while studying the cosmological 
evolution of scalar-tensor theories at earlier epochs in the matter-dominated era: The theory 
is attracted towards a maximum of F [i.e., a minimum of In A(ip)) so that the present value 
of |cto| is expected to be extremely small. Finally, this figure also displays the maximum 
value of |«o| for which the parametrization Z(&) = 1 of action ( |2.1| ) has a meaning. Beyond 
|a | — l/v3 (i.e., for a Brans-Dicke parameter — | < u < 0), one would get $q 2 < in this 



parametrization. In other words, Eqs. (|5.3j)-(|5T4|) cannot be integrated consistently beyond 



z ~ 1.58 if one sets Z($) = 1, whereas the Brans-Dicke or the Einstein-frame representations 
show that the theory can be mathematically consistent up to z ~ 3.5 (if' 2 remains positive). 
This underlines that the Z = 1 parametrization may be sometimes pathological. 

Our numerical integration of Eq. (|5.6b|) not only allowed us to check the positivity of 



the scalar field energy, but also to reconstruct parametrically the matter-scalar coupling 
function A(ip). Since A = F -1 / 2 , Eq. ( |2.4c|) , we know that A(z) is finite and strictly positive 



over the interval [0,2 max [, but we also checked that it is single valued over this interval. 
This means that if <p{z) can take several times the same value for different z, they must 
correspond also to the same value of A(z). Actually, since Eq. ( |5.6b| ) does not fix the sign 
of dcp/dz, one should keep in mind that tp can oscillate around a constant value <£> m i n . If 
the numerical integration confuses the two points (p m { n ± s, but if A(<p) happens not to be 
symmetrical around <p m m, it may look like a bi- valued function. When such a situation 
occurred in our programs, we always verified that a single- valued A((p) could be defined 
consistently by unfolding it around the oscillation points of tp. Figure 3 illustrates such a 
situation, for an intentionally unrealistic value of |ao| i n order to clarify the plots. [The 
value |ao| = 1 is inconsistent with the solar-system bound (|3.2j), but it corresponds anyway 
to a mathematically consistent theory, although the Z = 1 parametrization cannot be used 
in this case.] 

All the functions In A(ip) that we reconstructed have similar convex parabolic shapes. 
This is consistent with the results of Refs. |2"5] , |2l)[ , showing that the scalar field is generically 
attracted towards a minimum of In A(<p) during the expansion of the universe. If we had 
found models such that the present epoch (z = 0) is close to a maximum of In A, this would 
have meant that the theory is unstable, and that we have extremely fine tuned it to be 
consistent with solar-system constraints. On the contrary, the convex functions In A((p) that 
we obtained show that these scalar-tensor models are cosmologically stable, i.e., that the 
tight bounds (|3.2|) are in fact natural consequences of the attractor mechanism described in 
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We have checked that reducing the parameter f2A,o allows us to extend the integration 
region in the past, consistently with the above analytical results. For instance, when we 
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vary f2\.,o, still satisfying Q m ,o = 1 _ ^a,o and setting f2 Kj0 = 0, we find that fi\ t o < 0.02 is 
required in order to integrate the equations up to a redshift z = 5. This would correspond 
to fi m ,o/^A,o > 50, i.e., 100 times larger than present estimates. 

We also added random noise to our assumed H(z), Eq. ( |6.1|) , and verified that the 
conclusions are not changed qualitatively provided H(z) is known over a wide enough redshift 
interval. This means that the experimental determination of the luminosity distance Dl(z) 
needs not be very precise to be quite constraining, provided redshifts of order z ~ 2 are 
probed. As an illustration, let us take the exact expression H(z) of Eqs. (|6.1|) - (]6.2|) for 
discrete values of the redshift, say z — 0, 0.1, 0.2, 0.3,..., and let us add or subtract randomly 
between and 30% to the corresponding H(z). Then, we may fit a polynomial through 
these "noisy" values of H(z), and use our numerical programs to integrate the background 
equations (|5.3|) - ([^) and reconstruct F. We found that there always exists a maximum 
redshift beyond which F is negative (and the theory thus inconsistent). Figure 4 displays 
the two extreme values of z ma , x that we obtained with hundreds of such "deformed" H(z): 
It is sometimes even smaller than for the "exact" H(z) of Eq. fl6.1|) , and sometimes larger 
but never greater that ~ 2. It should be noted that for the 30% noise we chose, the H(z) 
of pure GR with a vanishing cosmological constant could have been obtained. In that case, 
a potential-free scalar-tensor model with $ = const, would of course have fitted perfectly 
this H(z) up to z — > oo. The reason why we never managed to go beyond z max ~ 2 is 
that we considered random noise, instead of such a precise bias of our assumed function 
H(z), Eq. ( |6.1| ). We are aware that our deformed functions of Fig. 4 do not reproduce a 
realistic experimental noise. However, they illustrate in a well defined way that an inaccurate 
determination of H(z) over a wide redshift interval is actually more constraining than a 
precise measurement over a small redshift interval only. 

The conclusion of the present subsection is therefore that a scalar-tensor theory without 
potential can accommodate a Hubble diagram consistent with ( |6.1| ), but only on a small 
redshift interval if £Ia,o is significant. The experimental determination of the luminosity 
distance Dl(z), either accurately for z ;$ 1 or even with large (tens of percents) uncertainties 
up to redshifts z ~ 2, severely constrains this subclass of theories. Future observations should 
thus be able to distinguish them from general relativity, and to confirm or rule them out 
without any ambiguity. 

It is worth noting that such future determinations of Dl(z) would a priori be much more 
constraining than solar-system experiments and binary pulsars tests. Indeed, although the 
precision of the latter is quite impressive (see e.g. |I^,[n],f2t],^I[] ) , they anyway probe only the 
first two derivatives of hi A(ip), Eqs. ( |3.2|) - (|3.3|) , whereas cosmological observations should 
give access to the full shape of this function. 

Let us also recall that the constraints we found crucially depend on the fact that the 
theory should contain only positive-energy excitations to be consistent, and notably that the 
function F should remain always strictly positive. We did not use any other cosmological 
observation, but obviously, once the microscopic Lagrangian of a scalar-tensor theory has 
been reconstructed using D L (z), all its other cosmological predictions should also be checked. 
For instance, a bound F nnc > 0.86 F is given in Ref. |34] for the value of the function F at 
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nucleosynthesis timeEl If one assumes that F(z) is monotonic, the reconstructed function of 
Fig. 1 would not be consistent with this nucleosynthesis bound beyond z ~ 0.3. This would 
be even more constraining than the bound z < 0.68 we obtained just from mathematical 
consistency requirements. Alternatively, a reconstructed function F(z) like the one of Fig. 1 
would be consistent with the above nucleosynthesis bound only if it were non-monotonic 
beyond z ^ 0.6. Although this would not be forbidden from a purely phenomenological 
point of view, this would be anyway unnatural, and more difficult to justify theoretically. 

B. Massless scalar field and an (arbitrary) nonzero cosmological constant 

To confirm the results of the previous subsection, let us now consider the case of a 
massless scalar field together with a cosmological constant whose value differs from the one 
entering our assumed H(z), Eq. (|6.1| )- (|6.2| ). The question that we wish to address is the 
following: Can part of the observed Qa,o be due to the presence of a massless scalar field? 

To impose a cosmological constant in a scalar-tensor theory, one would naively choose a 
constant potential U($>) in action ( |2.1| ). However, as shown by Eq. ( 2.4d| ), the correspond- 



ing potential V(p) of the helicity-0 degree of freedom ip would not be constant in this case 
(because -F($) is a priori varying), and its second derivative would give generically a non- 
vanishing scalar mass. To avoid any scalar self-interaction, and in particular to set its mass 
to 0, one needs in fact to impose V(<p) = const, in the Einstein-frame action fl2.5| ). This 
defines a consistent cosmological "constant" in a massless scalar-tensor theory. Note that 
the corresponding Jordan-frame potential U($>) is then proportional to F 2 ($), and therefore 
that it does not correspond to the usual notion of cosmological constant in action ( |2.1|) . 

Since our assumed "observed" H(z) involves a parameter denoted fi\,o, Eqs. fl6.1|) - flOD , 
let us introduce a different notation for the contribution due to the constant potential V : 

<V,o = W . (6.H) 

It is easily checked that for Q v ,o — ^a,o> the solution A(<p) = 1 (or F($) = 1) is recovered, 
i.e., a scalar field minimally coupled to gravity with a constant potential acting like a 
cosmological constant. Indeed, in terms of the function f(x) defined in ( |6.3Q , Eq. ( |5.3|) reads 

(tt Afi + n mfl x 3 )xf"(x) + ln mfl x 3 f(x)-4n Kfi xf(x)-6n Vfi xf 2 (x) = m mfi . (6.12) 



Note that this is now a non-linear equation in F, contrary to Eq. (|6.4j) above for the case of 
a vanishing potential. If Qv,o — ^a,o, one finds that f(x) = x~ 2 is an obvious solution, i.e., 
F(z) = F Q = const. A constant scalar field $ (or ip) then satisfies Eqs. (|5.4| ) -(|5.7| ) . 

If we now consider a scalar-tensor theory for which f2y,o differs from the "observed" 
Oao ~ 0.7, we find that like in the previous subsection, there exists a maximum redshift 



5 See however Ref. p7|] , in which extremely small values of F nuc /Fo = Ag/A^ uc are shown to be 
consistent with the observed abundances of light elements, provided d 2 A(ip) / dtp 2 is large enough, 
where A(<p) is the matter-scalar coupling function (2.4c). 
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^max beyond which F(z) becomes negative, and therefore beyond which the theory loses its 
mathematical consistency. Figure 5 displays this maximum redshift as a function of Qy,o- 
We plot this figure for the initial condition Fq = (i.e., for a theory which is presently 
indistinguishable from GR in the solar system), but as before, we verified that the curve 
is almost identical if one takes the maximum value of \Fq\ consistent with the solar-system 
bound (^3). We also assume f2 K)0 = (spatially flat universe) for this figure, as we know 
from the previous discussion that a value even as large as |fi K ,o| ~ 0.2 does not change 
qualitatively the results. 

For fl v ,o — 0, we recover the result z max « 0.66 derived above for a vanishing potential. 
When Qy,o < 0, -2 max becomes even smaller. As expected, this is worse than in the potential- 
free case. On the contrary, when flv,o is positive (i.e., when it contributes positively to part 
of the "observed" ^a,o)> the maximum redshift ;z max increases. This is just due to the fact 
that our massless scalar field needs to mimic a smaller fraction of the "observed" Q\ t o, so 
that the theory can remain consistent over a wider redshift interval. However, we find that 



z max is still smaller than 1.5 for Qv,o < 0.6, and a H(z) of the form ( |6.1| )- fl6T2j ) observed up 
to z ~ 2 would thus suffice to rule out the model. If such a H(z) could be confirmed up to 
z ~ 5, one would need Qv,o > 0.694 for our massless scalar-tensor theory to fit it! Even so, 
the theory would anyway become pathological at slightly higher redshifts. In conclusion, a 
massless scalar cannot account for a significant part of the observed cosmological constant 
if H(z) is experimentally found to be of the form (|6.1|) over a wide redshift interval. 

Let us note finally that for Qy,o > ^a,o, Eq. ( |6.12|) does admit strictly positive solutions 
for / (or F) up to arbitrarily large redshifts. This ensures that the graviton energy is 
always positive. However, it is now the scalar field which needs to carry negative energy. 
Indeed, Eq. (|5.7|) gives a negative value for ip' 2 , basically because of the presence of the large 
negative number —U in this equation. [In the Z{<&) = 1 parametrization, $ /2 is obviously 
also negative, because of Eq. ( |5.6a| ), or directly from Eq. ( |5.5|) which also involves a — U 
term.] 

Therefore, there is only one possibility for a consistent massless scalar-tensor theory to 
reproduce (|6.1| ) over a wide redshift interval: It must involve a cosmological constant, whose 
contribution Qv,o is equal to (or very slightly smaller than) the parameter Q\ t o entering fl6.1| ). 
In other words, the theory should be extremely close to GR plus a cosmological constant, 
and the massless scalar field must have a negligible contribution. This illustrates again the 
main conclusion of our paper: The experimental determination of the luminosity distance 
Dl{z) over a wide redshift interval, up to z ~ 2, will suffice to rule out (or confirm) the 
existence of a massless scalar partner to the graviton. 

C. Reconstruction of the potential U from a given F 

In the previous two subsections, the matter-scalar coupling function F($) (or A((p)) was 
reconstructed from the assumed knowledge of H(z), for theories whose potential U($>) (or 
V(<p)) had a given form. We now consider the inverse problem. We still assume that future 
observations will provide a Hubble diagram consistent with ( |6.1|) -( p\2"D , but we wish now to 
reconstruct the scalar-field potential U for given forms of the coupling function F. 
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1. Generic scalar-tensor theories 

We first consider a generic two-parameter family of scalar-tensor theories, which has 
already been studied in great detail for solar-system, binary-pulsar and gravity-wave experi- 



ments [pl| , p3[ , as well as for cosmology starting with the matter-dominated era |25| and even 
back to nucleosynthesis PT| . Its definition is simplified if we work in the Einstein frame 
( p.4|) - (p.5|) . The matter-scalar coupling function is simply given by 

In A(ip) = a (ip - <p ) + -po(<p - (f ) 2 , (6.13) 

in which the present value of the scalar field, ipo, may be chosen to vanish without loss of 
generality. Any analytical function In A(ip) may be expanded in such a way, but we here 
assume that no higher power of cp appears, i.e., that h\.A(<p) is strictly parabolic: It depends 
only on the two parameters «o and /3q. The latter is a simplified notation for (da/dip)o, and 
should not be confused with the post-Newtonian parameter (3 defined in ( [3.1b] ). [Actually, 
this equation shows that (3 ~ l + |«o/^o-] Solar-system experiments impose |«o| < 1-4 x 10~ 2 , 



Eq. (|3.2|) , while binary pulsars give (3q > —4.5, Eq. (|3.3|) , for this class of theories. We first 
study these models for the case of a spatially flat universe (fl^o = 0)- 

As shown by Eq. (|4.7|) , a constant scalar field ip = (p may be a solution if a = 
(so that a((p) ex (<p — (p ) vanishes too) and if the potential V((p) is also constant. Our 
assumed H(z), Eqs. ( |6.1|) - (|6~2] ), can thus always be reproduced if the parameter a vanishes 
identically, and the reconstructed potential merely reduces to the constant V = ~HqQa,o- 
This corresponds simply to GR plus a cosmological constant, and the massless scalar degree 
of freedom ip remains unexcited, frozen at an extremum of the parabola ( |6.13j ). Actually, 
Eq. fl4.7| ) shows that this extremum corresponds to a stable situation only if it is a minimum, 



i.e., if (3q > in (|6.13|) . This is consistent with the results of Refs. p5| , p6| : If the theory 



involves a cosmological constant whose value equals the "observed" one in Eqs. ( |6.1| )-( |BT2] ), a 
massless scalar field is cosmologically attracted towards a minimum of the coupling function 
\a.A{ip), and the present value of its slope, «o, is expected to be generically very small. 

On the other hand, if «o is not assumed to vanish, say if its value is comparable to the 
solar-system bound (|3.2| ), then our reconstruction of the potential V(ip) from Eqs. (|5.3|) - 
( p.6| ) leads to serious difficulties. Their nature depends on the magnitude of the curvature 
parameter /3 of parabola ( |6.13| ). 

String-inspired models [^TJ suggest that j3 may be as large as 10, or even 40. With such 
large values (and assuming non- vanishing a ), our numerical integrations of Eqs. (|5.3|) - (|5.6|) 
give concave potentials V(<p), unbounded from below. This corresponds to unstable theories, 
and thereby to extremely fine-tuned initial conditions: Changing slightly the derivative of 
the scalar field, dip/dz, at high redshifts would a priori yield a totally different universe 
at present. This result tells us that this kind of models cannot be consistent over a wide 
redshift interval with the exact form of H(z) we chose in (|6.1| ), unless the parameter a® is 
extremely small. Actually, this is just another way to present the results of Refs. PREBJ: 
Since they predict that «o should be almost vanishing at present, assuming a significant 
non-zero value implies that the theory is unnatural. 

To obtain convex-shaped potentials V(tp) (i.e., stable theories) while still assuming a non- 
vanishing a , we typically need values of |/3 | ~ 4. However, the reconstructed potentials 
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always exhibit sudden changes of their slope. Basically, they reproduce a cosmological 
constant over a finite interval around (fio (i.e., around z = 0), and become rapidly divergent 
beyond a critical value of the scalar field (depending on /3 ). Therefore, as in subsection 
VLB above, we find that such scalar-tensor models can reproduce ( |6.1|) only if they involve 
a cosmological constant, whose energy contribution is close to the parameter Qa,o entering 
H(z), and if the scalar field has a negligible enough influence. In other words, such models 
would not explain the small but nonzero value of the observed cosmological constant by 
a "quintessence" mechanism, and would not be more natural than merely assuming the 
existence of A. 

The above results are significantly changed if we take into account the possible spatial 
curvature of the universe. Indeed, smoother potentials V(ip) are obtained for closed universes 
(fl K fi < 0), and the present value of the cosmological constant thus becomes more "natural". 

To illustrate this feature, let us consider the case of a minimally coupled scalar field (as 
in fl4]|), corresponding to a = /3 = in Eq. ( |6.13| ). For an open universe (fi K ,o > 0), we 
find from Eq. ( |5.6| ) that the scalar field would need to carry negative energy to reproduce 
d6.1| ). On the other hand, for a closed universe (Q K) o < 0), one can derive analytically the 
parametric form of the potential V((p). It can be expressed in terms of the hypergeometric 
function 2 ^i (a, b] c; x) (solution of the differential equation x(l — x)T" + [c— (a + 6+l)x]JF' — 
abF = 0): 

V= (ln Afi + n Kfi x 2 )H 2 , (6.14a) 



,2 

^ = ±x F^^ 1 l3'2 ; 3 ; -a^ x J ' (6 - 14b) 

where as before x = 1 + z. If |f2 Kj0 | is very small, we recover that V((p) exhibits a sudden 
change of slope, as was obtained above in the flat case. This is illustrated by the left panel of 
Fig. 6. On the contrary, if |f2 K0 | is large enough, the same analytical expression ( p. 14 ) gives 



nice regular potentials, like the one displayed in the right panel of Fig. 6. This reconstructed 
V(<p), as well as those obtained numerically for weakly varying In. A(<p), Eq. ( |6.13|) , are 
natural in the sense that they can be approximated by the exponential of simple polynomials 
in if. In that case, the observed value of the cosmological constant does not appear as a 
mere parameter introduced by hand in the Lagrangian, but corresponds basically to the 
present value of 2V((po). It should be noted that a value as large as f2 K)0 = —0.1 is not 
excluded by the latest Boomerang data, though it would be problematic in the framework 
of the inflationary paradigm. 

In conclusion, the existence of non-singular solutions over a long period of time is again 
the constraining input. A non-minimally coupled scalar field is essentially incompatible with 
( |6.1| ) over a wide redshift interval, unless the scalar field is frozen at a minimum of In A(ip) 
(consistently with |25|,|2B|]). If future experiments provide a Hubble diagram in accordance 
with ( |6 . 1| ) and also give a very small value for Q Kj o, it will be possible to conclude that scalar- 
tensor theories (either non-minimally or minimally coupled) cannot explain in a natural way 
the existence of a cosmological constant. On the other hand, if the universe is closed and 
\ft Kj o\ large enough, a "quintessence" mechanism in a scalar-tensor theory seems more natural 
than a mere cosmological constant. 
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2. Scaling solutions 

The above conclusions can be confirmed by starting from a given F(z) (or A(z)), rather 
than -F($) (or A(ip)). We consider here "scaling solutions", i.e., we assume that these 
functions behave as some power of the scale factor a. One may for instance write F(z) = 
(a/a ) p = (1 + z)~ p , with p > 0. As before, our aim is to reconstruct a regular potential 
V(<f>) from the knowledge of H(z), assumed to be of the form (p.l|)-(|6T2D. 

The strongest constraint on this class of theories is imposed by the solar-system bound 
(p.2|). Indeed, using the definition ( |2.7| ) for a(<p), one can also write it as a(<p) = —F'/(2ip'F), 
and Eq. ( |5.6|) evaluated at z = then yields the following second-order equation for p : 

(1 - a 2 )p 2 - (2 + 3n m , )a 2 p + A^ a 2 = . (6.15) 



Note that this equation does not depend on the full form of Eq. ( pM| ), but only on its first 
derivative at z — 0, i.e., on the deceleration parameter q = (H'/H) — 1. The constraints 
on p derived below are thus valid as soon as go is of order ~ — |, consistently with the 
estimated value (6.2) for fi m ,o- 

In the case of a spatially flat universe (fi K) o = 0), Eq. ( |6.15| ) gives immediately p = 
(2 + 3Q mi o)al/ (1 — al) ~ 3al, so that the solar-system bound (|3]2|) imposes p < 6 x 10~ 4 . 
Therefore, the scalar field needs to be almost minimally coupled. If p vanishes identically, we 
recover as before the trivial solution of GR plus a cosmological constant, together with an 
unexcited minimally-coupled scalar field. On the other hand, if p does not vanish, one finds 
that the scalar field needs to carry negative energy beyond z ~ 1.4. Even without trying 
to reconstruct the potential V(ip), one can thus conclude that such scaling solutions would 
be ruled out by the observation of a H(z) of the form (|6.1|) up to z ~ 2. Paradoxically, 
this result is valid even for an infinitesimal (but nonzero) value of p. Indeed, there exists 
a discontinuity between the case of a strictly constant F and that of a scaling solution 
F(z) = (1 + z)~ v . At first order in p, and still assuming f2 K)0 = 0, one can write Eq. ( |5.6|) as 

\2 12 _ o (\ ft i N 1\ ^m,o(l "~ Z Y 



2(l + Z y^=p-3p(Hl + z)-- ^ /, ,3 +0^). (6-16) 

This equation confirms that if' 2 — > when p — > 0, and therefore that the scalar field tends 
towards a constant in this limit. However, it carries positive energy (<// 2 > 0) only if 

6/ 3f2 mj o 



1 + Z y (\n(l + z)-- )<— ^L. (6.17) 



Since the right-hand side is estimated to be ^5 1, the maximum value of z is obtained for 
ln(l +z) w 5/6, so that the large numerical factor coming from (1 + z) 3 in the left-hand side 
is compensated by the small term inside the second parentheses. Working iteratively, this 
maximum redshift can be better approximated by z max ~ e 5//6 — l + (f2A,o/3fi mi o)e~ 5//3 ~ 1.45, 
and the actual numerical resolution of equality ( |6.17| ) for the values ( |6.2| ) gives z max = 1.429. 



Therefore, even if p is vanishingly small, a scaling solution F(z) = (1 + z)~ p cannot be 
consistent with ( |6.1| ) beyond this maximum redshift. This illustrates once more that the 
experimental determination of H(z) up to z ~ 2 would be more constraining that solar- 
system experiments for this class of theories, provided one takes into account the requirement 
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of positive energy. Let us underline that the above value for z max is valid for a monomial 
F = (a/a ) p but not for more complicated polynomial expressions. Indeed, as shown for 
instance in Sec. VLB above, there do exist scalar-tensor theories consistent with ( |6.1| ) up 
to arbitrarily large redshifts, and they do not need to be strictly equivalent to GR plus 
a cosmological constant (although they must be close enough to it). Moreover, the above 
maximum redshift is a consequence of the exact form for H(z) we chose in Eq. ( |6.1|) . A 
slightly different function may of course allow a positive-energy scalar field up to much 
higher redshifts. It suffices that the right-hand side of Eq. ( |5.6b|) be strictly positive for 
F ~ const., and the case of a closed universe discussed below provides an example, since 
the contribution — Q, K)0 is then positive in Eq. ( |5.6b| ). 

The case of a spatially open universe (f2 Kj0 > 0) is forbidden by Eq. (|6.15|) , unless Q K ^ is 



smaller than ~ ^o < 10 . Such a situation would be indistinguishable from the spatially 
flat case. 

In a spatially closed universe (£l Kt o < 0), p is given by the positive root of the second- 
order equation Eq. ( |6.15|) . Remembering the solar-system bound (|3.2| ), one may consider 



the case «q -C |^ re ,o|, and one gets p ~ 2\a \J— Q K fl- Even if one considered values of Q K} o 
as large as —0.1, this would limit p to ~ 10~ 2 . Therefore, in this case again, solar-system 
constraints impose that the scalar field should be almost minimally coupled, if one looks 
for such scaling solutions. The difference with the spatially flat case is that Eqs. (|5.3|)-(Jx7D 
can now be integrated for any redshift z (from future infinity, z = —1, to arbitrarily large 
z). Since F(z) needs to be almost constant, we recover solution ( |6.14| ) for the potential 



V((p). As in Sec. VI.C.l above, we can thus conclude that such models would be consistent 
with ( |6.1| ) over a wide redshift interval only if they are (almost) minimally coupled, and 
they would provide a natural "quintessence" mechanism to explain the presently observed 
cosmological constant only if the universe is (marginally) closed. 

Let us end this paragraph by a remark concerning scaling solutions, for which the scalar- 
field energy density scales like a power of a. As mentioned in the Introduction, they have 
attracted a lot of attention recently. For a minimally coupled field, the possible scaling 
behaviors and the corresponding potentials can be classified |L2|. As for a non- minimally 
coupled field, a subclass of theories was considered in pE5|| , for which 

[/($) = OF($) M , (6.18) 

where C and M are constants. Since, besides these two constants, there is only one unknown 
function of $, the knowledge of H (z) suffices to reconstruct the full microscopic Lagrangian 
from Eqs. (|5.3|)-(|5!6|) above. However, the main conclusion of Ref. p5[ can be recovered 



from a simple argument, without any numerical integration. Indeed, it was shown in this 
reference that there exists a universal behavior of these theories, depending on M but not 



on the precise shape of F(&). As emphasized in [f45|, this result was obtained in the strong 
coupling limit, corresponding formally to Z(§) — > in action ( |2.1| ). Taking into account the 
assumed relation (|6.18| ), the class of theories under consideration is thus defined by 

S = j^- J d 4 xV=g(Fm R - 2CT($) M ) + S m [ij mi g, v ] . (6.19) 

If we now introduce a new scalar variable \1/ = F(<&), we notice that $ disappears totally 
from the action. No physical result can thus depend on the precise form of F(<&), and we 
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recover the conclusion of j|5| . The constant C may also be set to 1 by a change of length 
units, and this class of theories is thus parametrized by the single real number M. Any 
physical prediction must therefore depend only on M. 

VII. CONCLUSION 

In this work we have investigated the constraints that arise from the experimental knowl- 
edge of the luminosity distance in function of the redshift up to z ~ 2, corresponding to 
H(z) given by (|6.1| )- (|6.2| ). In particular, our universe is then presently accelerating and we 
have studied the viability of subclasses of scalar-tensor theories of gravity. We have shown 
that the subclass of models in which the scalar partner $ of the graviton has no potential 
at all, and which satisfy the present-day existing constraints, are inevitably ruled out if an 
expansion of the form of Eq. (|6.1|) holds even for a redshift interval as tiny as z < 2 (see 
the precise numbers in section VI). We see that these theories become pathological in the 
form of a vanishing F, already at such low redshifts for which H(z) will be experimentally 
accessible in the near future [|TjJ . Hence we show that a cosmological observation of the back- 



ground evolution according to Eq. (|Q|) in the "recent" epoch will be enough to rule out such 
models. [On the other hand, future observations might provide a H(z) which confirms the 
existence of a scalar partner to the graviton and rule out pure GR!] The main reason why we 
obtained so constraining results is that we took into account the mathematical consistency 
of the theory, i.e., the fact that it should contain only positive-energy excitations to be well 
behaved. This requirement severely restricts the class of viable models. 

A non-flat universe can alleviate in some cases the tight constraints we found. However, 
the latest CMB data released by Boomerang and Maxima ||^]] favor a flat universe (in 
accordance with the inflationary paradigm), and only a marginally closed universe is still 
allowed by the location of the first acoustic (Doppler) peak at I ~ 200, while an open 
universe is more unlikely. 

The most impressive conclusion is that future cosmological observations may prove to 
be more constraining for massless scalar-tensor theories than solar-system and binary-pulsar 
tests. Indeed, even if the determination of the luminosity distance Dl(z) will not reach very 
quickly the impressive accuracy obtained in the solar system or with binary pulsars, it will 
nevertheless give access to the full coupling function -F($) in action ( |2.1| ), or A(cp) in the 
Einstein-frame rewriting (|2.5|) , whereas only its first two derivatives are presently probed. 
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FIGURES 
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FIG. 1. Reconstructed F(z) [i.e., Brans-Dicke scalar $bd(^)] and Einstein-frame scalar (p as 
functions of the Jordan- frame (i.e., observed) redshift z, for the maximum value of |qo| allowed 
by solar-system experiments, and for a vanishing potential. The helicity-0 degree of freedom (p 
diverges at z max ~ 0.68. 
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FIG. 2. Maximum redshift z consistent with the positivity of energy of both the graviton and 
the scalar field, as a function of the parameter |«o|- This figure corresponds to the case of a 
vanishing scalar-field potential, and we fit the exact H(z) predicted by general relativity plus a 
cosmological constant (GR +A). 
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FIG. 3. Two versions of the reconstructed coupling function InA(ip) for |ao| = 1) the dashed 
one looking bi- valued, but the (single- valued) solid one giving the same predicted H(z). This 
figure still corresponds to the case of a vanishing scalar-field potential, and we fit the exact H{z) 
predicted by GR +A. 
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FIG. 4. Random deformations of the H(z) predicted by GR +A (with Q\fl = 0.7), and cor- 
responding maximum value of the redshift z consistent with the positivity of energy. The dashed 
lines indicate the region in which random points have been chosen at regular intervals of z. The 
thin solid lines correspond to two polynomial fits of such random points. Note that they can differ 
from the GR +A curve even more than the dashed lines. The dotted line labeled simply "GR" 
corresponds to a vanishing cosmological constant A. Such a bias of the GR +A curve changes z max 
much more that the random noise we considered. 
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FIG. 5. Maximum redshift z consistent with the positivity of energy, as a function of the value 
of a constant potential V (case of a massless helicity-0 degree of freedom </?). 
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FIG. 6. Minimally coupled model F = 1 in a spatially closed FRW universe, respectively for 
&k,o = — 10~ 3 (left panel) and Kj o = —0.1 (right panel). In both cases, the potential V(tp) is 
analytically given by Eq. (|6.14 ). Note that the reconstructed potential does not have a "natural" 



shape if |^ re ,o| is too small: The present value of Oa,o is n °t explained by a quintessence mechanism, 
and the corresponding scalar-tensor theory is basically equivalent to GR +A. On the contrary, if 
|Oa,o| is large enough, the potential has a nice smooth shape, and its present value (at ip — ipo = 
on the figure) basically corresponds to the observed Oa,o- 
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